A Dynamical Study of the Disordered Quantum p = 2 Spherical Model 
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We present a dynamical study of the disordered quantum p — 2 spherical model at long times. Its 
phase behavior as a function of spin-bath coupling, strength of quantum fluctuations and temper- 
ature is characterized, and we identify different paramagnetic and coarsened regions. A quantum 
critical point at zero temperature in the limit of vanishing dissipation is also found. Furthermore 
we show analytically that the fluctuation-dissipation theorem is obeyed in the stationary regime. 



I. INTRODUCTION 

The dynamics of quantum systems far from equilib- 
rium present many conceptual challenges, particularly in 
the area of quantum glasses. The classical analogues 
of these systems display broad relaxation spectra; fur- 
thermore their long-time behavior depends explicitly on 
details of sample history. In this paper a motivating 
question is whether the introduction of quantum fluctu- 
ations in a classical system out of equilibrium results in 
qualitatively new behavior at long time-scales. To date, 
the introduction of quantum fluctuations in theoretical 
glassy models has not led to qualitative changes, e.g., 
new slow quantum modes, in their long-time dynamics 
IH EH3| Quantum fluctuations often drive the glass 
transition to be first-order 

0111111113, 

and thus only 

qualitatively affect the dynamical behavior on short-time 
scales. However there are examples where, in the ab- 
sence of dissipation, quantum fluctuations lead to a quan- 
tum critical point at zero temperature; examples include 
the periodic Josephson array (Tj and infinite-component 
quantum rotors 0, ^| . Detailed characterization of the 
phase space as a function of dissipation, temperature, 
and quantum fluctuations remains to be studied in both 
cases. 

Careful experimental studies of the dipolar magnet 
LiHo x Yi_ x F4 indicate that the current theoretical pic- 
ture of quantum glassy systems is not the full story |l4j . 
An appropriate model for the description of the mag- 
netic properties of this material consists of Ising spins, 
representing the states of the magnetic ions, residing on 
random sites, and interacting via short-range ferromag- 
netic couplings and long-range dipole-dipole interactions. 
The sign of the latter interaction depends on the relative 
position of the spins, leading to the possibility of com- 
peting interactions and to frustration. For x = 1, the 
short-range component of the interaction dominates and 
the system is a ferromagnet. However, the Curie temper- 
ature decreases as x is reduced. The effects of quantum 
mechanics can be tuned in this system by application 
of a magnetic field transverse to the spins. Even in the 
disordered ferromagnetic state (a; ~ 0.4), there are qual- 



* Present address: Department of Physics, Rutgers University, Pis- 
cataway, NJ 08855, USA 



itative effects of the transverse field, particularly near 
the observed quantum critical point, that have not been 
described by theoretical models [lj, [l5|, . 

Motivated by these experiments, we study the dynam- 
ics of a long-range quantum disordered ferromagnet, the 
quantum spherical model with Gaussian random two- 
spin interactions (p = 2). Classically the statics of this 
model are trivial, as the free energy has only two min- 
ima, and it can be solved exactly with a replica sym- 
metric ansatz |l7j| . However, for almost any initial con- 
ditions, the dynamics of this classical system are out-of 
equilibrium with pow er-law decays in many of the ob- 
served quantities [lg ■ The origin of this slow relaxation 
is very different from that of classical mean-field spin 
glasses. In these latter systems the time evolution is de- 
termined by the complex structure of the phase space 
|l9| . By contrast the free energy landscape of the p = 2 
case is very simple. However, statistically the overlap 
between the initial and the equilibrium configurations is 
very small. This is qualitatively similar to the situation 
in a weakly disordered ferromagnet after a rapid temper- 
ature quench where the resulting domain structure bears 
little resemblance to the uniform equilibrium state. For 
this reason we refer to the non-ergodic phases of the p = 2 
spherical model as coarsened. 

Naturally the infinite-range nature of this Hamilto- 
nian implies the absence of any length-scale in real-space. 
The spherical character of the spins combined with the 
quadratic interaction lead to reduced nonlinearities in 
this disordered p = 2 model, even compared to other 
mean-field spin glasses. These two key features, specifi- 
cally the infinite-range coupling and the minimal nonlin- 
earities, make this model amenable to analytic treatment, 
and it can serve as a testing ground for necessary assump- 
tions made in more complicated glassy systems. The dis- 
ordered quantum p = 2 spherical model bears a strong 
resemblance to that of the infinite-range M-component 
quantum rotors in the limit M — > oo [T3 . Il3|| . More ex- 
plicitly, the M-component quantum rotor Hamiltonian is 
composed of M replicas which each individually resemble 
the p = 2 system, where the inter-replica coupling scales 
inversely with M [T^ |. Here we extend previous studies 
of these and related rotor Hamiltonians d IT2I Il3l| by 
performing a detailed characterization of the phase be- 
havior of the disordered quantum p = 2 spherical model, 
noting that our results recover those previously obtained 
in particular well-defined limits. 
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In this study we characterize the phase behavior of 
the disordered p — 2 spherical model as a function of 
quantum and thermal fluctuations and the spin-bath cou- 
pling; we identify its distinct paramagnetic and coars- 
ened regions by studying its long-time dynamical behav- 
ior. Known results are reproduced in both the high- and 
the zero- temperature limits. We also demonstrate ana- 
lytically that the fluctuation-dissipation theorem remains 
valid in the stationary regime. 

The plan of the paper is as follows. In Section [Q] 
we present the model. Its dynamical equations, which 
we derive using the Keldysh formalism, are discussed in 
Section IIIII The stationary parts of the retarded and 
the statistical Green functions, that we find in Section 
IIVI are then shown analytically to obey the fluctuation- 
dissipation theorem. In Section an expression for the 
Edwards- Anderson order parameter is presented with the 
resulting emergence of a quantum critical point. Details 
of the phase diagram are discussed in Section IVII for all 
the coarsened and the paramagnetic regions. We end 
with a discussion, Section lVllI where, after summarizing 
the main features of this work, we mention several open 
questions that emerge from our study. 



II. FORMULATION OF THE PROBLEM 

In this paper we study the dynamics of the disordered 
quantum p — 2 spherical spin glass model coupled to an 
Ohmic bath. There are N spins that can be written in 
vector form as s = (si, S2> • • • > sjv), and we will assume 
that N — > oo. The Hamiltonian of the system takes the 
form 



H = Ho + Tiss + Hb + Hsb, 



(1) 



where Ho, Hss, Hb and Hsb are the free, the spin-spin, 
the bath, and the spin-bath contributions respectively. 

In the absence of interactions the free Hamiltonian of 
the spin system is given by 



Ho 



7T 

2m' 



(2) 



where the conjugated momenta ir = (wi, 7T2, . . . , ttn) 
obey the usual commutation relations [si, Sj] = [ivi, ttj] = 
0, and [7Tj,Sj] = —iSij; here we use a unit system where 
h = kg = 1 . The "mass" m is a constant that introduces 
dynamics into the problem. As we will explain below, the 
dynamics of the system can be tuned from the classical 
to the quantum regimes by varying the mass m. We note 
that there are other methods of quantization (e.g. appli- 
cation of a transverse field) besides this "kinetic" one' if 
done correctly, they should all yield the same results [6| • 
The spins interact with each other via a spin-spin in- 
teraction which, for the p = 2 model, is given by 



H SS [J] = 



N 
il<i2 



(3) 



with randomly distributed matrix elements Jy , that obey 
J,, = 0, JsaJuu • / , • /, = 0, and J 2 = J 2 IN, where ~ 
represents averaging over the random spin-spin interac- 
tion. 

In the spherical model the spins, Sj, are constrained 



to a spherical surface such that 



N 



»2\ 



Y . where 



(• • •) represent quantum averaging. Thus Sj are contin- 
uous variables, that behave like coordinates. However, 
since this model has been introduced as an analytically 
accessible spin problem, we will keep the notation and 
the nomenclature. The spherical constraint is initially 
introduced into the Lagrangian, which results in the ad- 
ditional term, ^z(t)s 2 , in the Hamiltonian. The time- 
dependent function, z(t), is a Lagrange multiplier. Thus 
the free Hamiltonian Ho (Q), is effectively replaced by 



Ho 



TV 

2m 



1 



z{t) s 2 



(4) 



Each spin in the system is connected to a system of 
Nb independent harmonic oscillators that constitute the 
bath. The free Hamiltonian of the bath is given by 



Hb 



It, 



(5) 



where x a = (x a i, x a 2, ■ ■ ■ , Xo,n) are the oscillator coordi- 
nates of mode a, and x a i is the coordinate of the oscilla- 
tor that interacts with spin s,. The conjugate momenta 
Pa = (P a i,P a 2, ••••Paiv) are related to x a through the 



commutation relations 



] = [P ai ,P aj ] = 0, and 



-iSij. The mass and frequency of mode a 



of the oscillators are given by M a and u) a , respectively. 
For simplicity we assume a linear interaction between the 
spins and the bath [2(| 



Hsb [Co] = ^2 CaXa ' S ' 



(0) 



where the coupling coefficients are given by C a . We are 
aware that this model for a bath can lead to unwanted 
divergencies due to low-frequency contributions in the 
absence of an infrared cut-off |2l| . However, this problem 
does not arise in our case due to the structure of the 
relevant spectral density, as we discuss later in the paper. 
Such a simple interaction between the spin system and 
the bath allows us to solve the problem analytically. 

The Hamiltonian for the spin system is actually that of 
a harmonic oscillator. Because the coupling to the bath 
© is linear in the spin operator, it does not alter this 
statement. In the paramagnetic phase the spin-spin in- 
teraction is negligible, and we recognize that y/z(t)/ 771 IS 
the eigenenergy of the oscillator. We are interested in the 
long-time behavior of z(t) when it saturates to a constant 
Zoo ■ In order to understand whether the dynamics of the 
spin system are governed by classical or quantum fluc- 
tuations, we compare the important energy-scale of the 
system to the temperature T . At constant T when m 
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is increased, the eigenenergy of the spin Hamiltonian de- 
creases. Once the dominant energy-scale is smaller than 
T, the system dynamics are governed by classical behav- 
ior. Similarly, when m is decreased, T becomes smaller 
than the dominant energy-scale, and quantum dynamics 
become more important. As a cautionary aside, we note 
that Zoo may depend on m in the paramagnetic regime; 
however we will show that to leading order it is either 
independent of m or proportional to m _1 , so that the 
preceding heuristic argument is still valid. 

In the coarsened phase the spin-spin coupling cannot 
be neglected. After averaging over this interaction, it 
is harder to recognize the eigenenergies. However, the 
results presented here indicate that when the coupling 
to the bath is infinitesimal, the signatory parameter is 
x = y mT 2 / J : for x large, the dynamics are dominated 
by classical effects whereas for x small quantum effects 
are involved. It is only at T — that the dynamics are 
purely quantum in nature. We conclude therefore that 
the energy-scales T and y/J/m should be compared in 
order to find whether the dynamics are purely classical 
or partially quantum (except at T = 0). This is con- 
sistent with the previous discussion of the paramagnetic 
phase, where small mass corresponds to a regime where 
quantum fluctuations could become dominant. 



III. THE DYSON EQUATIONS FOR THE SPIN 
GREEN FUNCTIONS 



The dynamics of the spin system can be characterized 
by Dyson equations associated with its Green functions. 
Here we use a nonequilibrium diagrammatic description, 
the Keldysh technique [23, |23| , and take the initial con- 
ditions to be random. Two Green functions, or rather 
three real functions, are needed to describe the system, 
they are the retarded and the statistical (or Keldysh) 
Green functions, that are defined by 



G T (t, = -iQ(t - t')( Si (t)si(t') - Si (t')si{t)), (7) 



and 



G S (t,t') = -i(8i(t)8i(f)+8i(tf)8i(t)), 



(8) 



respectively. For the p = 2 spherical model the Green 
functions are diagonal in the spin index, due to both the 
disorder averaging, and the infinite-range nature of the 
interaction. The retarded Green function is proportional 
to the "response" in the spin glass nomenclature, while 
the statistical Green function is proportional to the "cor- 
relation" . At equal times we demand that G s (t, t) = —2i 
due to the spherical constraint, thus setting z(t). 

When a system is at equilibrium the Green functions 
that describe it are a function of the time difference only. 
In that case one can Fourier transform G r (t — t') and 
G s (t — t') over the time difference to obtain G r {u>) and 
G s (lo). Furthermore, at equilibrium the retarded and the 



statistical Green functions are related via the fluctuation- 
dissipation theorem (FDT) 



(9) 



G s M = 2 lC oth^— JImG r M, 



where T is the equilibrium temperature of the system. 
When the spin system is in the paramagnetic phase, it is 
in equilibrium, and therefore it obeys FDT. Although in 
the coarsened phase the system is not in equilibrium, we 
will show that one can identify unique stationary contri- 
butions to the retarded and the statistical Green func- 
tions, that are related via the FDT, where the tempera- 
ture is that of the bath. 

When the spin system is in the coarsened phase, the 
Green functions depend on both times, and not only on 
the time difference. The spin Green functions in the time 
domain obey the Dyson equations 



d 2 



(10) 



6(t-t')+ / dt"Z r (t,t")G r (t",t'), 



and 



d 2 



G S (M') = (11) 

J dt" [£ s (i, t")G a (t", t') + E r (i, t")G s (t", t')} , 

where G a (t,t') = [G r (t',t)]* is the advanced Green func- 
tion, and "E T / S (t,t') are the self-energy functions. Note 
that the Dyson equations above apply to the system in 
the paramagnetic phase also, where z{t) is a constant 
and all other functions depend only on the time differ- 
ence. Assuming N is large, we consider the self-energy 
functions in the self-consistent-Born approximation [see 
Fig. Ufa)], in which case they are given by 

£ r / s (£, t') = D*' a {t - t') + J 2 G r/s {t, t'). (12) 

The bath contributions to the self-energy, which include 
the spin-bath interaction vertices, are D T / S (t — t'). We 
have neglected "crossing" diagrams [e.g., see Fig. ^b)], 
since they are at least a factor of 1/N smaller than the 
diagrams considered. This approximation is exact in the 
N — > oo limit. 

We remark that the equation for the retarded Green 
function (fTTJfl is completely decoupled from the statistical 
Green function. This is a result of two model-dependant 
features: (i) that p = 2, (ii) that the coupling to the bath 
is linear in the spin operators. In a p > 3 model, or in 
any p model with a nonlinear coupling to the bath, the 
equation for the retarded Green function would include 
the statistical Green function as well. 

The detailed nature of the relaxation of a spin system 
far from equilibrium depends on the time it has spent in 
its low-temperature phase before the start of the mea- 
surement. This phenomenon is known as aging. This 
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(a) 



(b) 
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G° 



G° D G 



G° G G 



FIG. 1: (a) Dyson equations depicted diagrammatically in a 
schematic manner, in the self-consistent-Born approximation. 
The dashed lines represent averaging over spin-spin interac- 
tion, (b) An example of a diagram with crossed spin-spin 
interaction averaging lines, that is proportional to 1/7V, and 
therefore is neglected in our calculation. 



dynamical behavior has been described using scaling and 
other phenomenological approaches [H^. |24| . More re- 
cently aging has been studied analytically in microscopic 
models where mean-field models are exact 0,01- ^ was 
realized that in these infinite-range systems the time evo- 
lution is dominated by motion in flat directions in phase 
space rather than by transitions over barriers; with ran- 
dom initial conditions, the system continues to evolve and 
never approaches its local equilibrium. For this "weak- 
ergodicity" breaking scenario [Hi. l24l I25ll26l| . it is crucial 
that the dimensionality of phase space is very large. Be- 
cause the spin system is out of equilibrium, its memory of 
its initial states will depend specifically on the time when 
the measurement commenced, t', and the time when it 
ended, t = t' + At, and not simply the time difference, 
At = t — t' . In the weak-ergodicity breaking scenario, 
if the waiting-time t' is finite, the system can evolve in 
the long-time limit so as to lose complete memory of its 
initial configuration at t' . This assumption of weak long- 
term memory has been shown to be self-consistent for 
several classical long-range glassy and coarsened models 
124 l25l l26j. In the quantum case, it is expected to 
hold in general 0, Q , though it may be violated near a 
quantum critical point where the decohcrcnce time di- 
verges. 

We assume that the weak-ergodicity breaking scenario 
applies to our system, and look for a solution that obeys 
it. More specifically, according to the weak-ergodicity 
breaking scenario, one expects the statistical Green func- 
tion to decay in the following manner 



lim^G s (M') = -2iq + G s ST (t-t'), (13) 
0, (14) 
0, (15) 



lim G| T (t-t') 



lim G s (t,t') 

t — >oo 



for t > t' . The Edwards- Anderson parameter, 
q = — lim t _ t /^oo limt/^oo G s (t, t')/2i, characterizes the 
coarsened phase. For large waiting times t' , the statisti- 
cal Green function is given by a term that is proportional 
to q, plus a stationary function G| T (i — t 1 ) that decays to 
zero at infinity. This serves as a definition of G| T (t — £'). 



Using the definition for G| T (t — f), the statistical 
Green function for any times t' and t can be written as a 
sum of G| T (i — t') and an aging function G s AG (t, t'), that 
is not time translationally invariant, 

G s (M') = Gl G (i,i') + G s ST (i-i'). (16) 

At equal times G s (t,t) — —2i due to the spherical con- 
straint, so that Eq. (|13|) at equal times yields 

lim G s AG (t, t) = -2iq and G| T (0) = -2i(l - q). (17) 

t — >oo 

In order to write the Dyson equations l|l() [l -l|ll jl . we 
need the contribution of the bath to the self-energy [see 
(|12p]. In the frequency domain the free-bath Green func- 
tions for mode a, that are associated with the free-bath 
Hamiltonian TLb, are given by 



w — u> a + i5/2 uj + uj a + i(5/2 ' 



(18) 



D s a (u) = -2ni coth (|M [S(uj - uj a ) + 5{lo + u a )] , 

where T is the bath temperature. Fourier transforming 
into the time domain one obtains 



D T a (t) = -2e(i)suW, 
D s a (t) = -2i coth cos to a 



(20) 
(21) 



We assume that the bath is in equilibrium, namely that 
it is so large that it is unaffected by its interaction with 
the spin system. Therefore we take the bath Green func- 
tions to be those in (|18l) and (|19fl . that depend only on 
the time difference. Thus the contributions of the bath 
to the self-energy are also time-translationally invariant. 
Furthermore, they are linear in the bath Green functions. 
They are obtained by multiplying the Green function of 
each mode by its density of states, (2M a u> a )~ 1 , and by 
the spin-bath coupling constant squared G\, and sum- 
ming over all modes 

/oo 
dwI(o;)sincjt, (22) 
-00 

for the retarded self-energy, and 

D s (t) = -2ij duj coth cos (ut), (23) 



for the statistical self-energy, where [2(j 

N b C 2 

^ 2M a w a 

a—l 

The bath temperature T enters the equations for the spin 
Green functions only through the bath contribution to 
the statistical self-energy. 



5 



In the case of an Ohmic bath I(oj) = 
(jui/tt) exp (— |w|/fi) [2jj, where 7 is a dimension- 
less coupling constant, and fl is an ultraviolet cutoff 
that will be taken to infinity later on. Thus for the case 
of an Ohmic bath 

^ = >^iW (25) 

When fl — > 00 one can Fourier transform l|25|) and obtain 
D T (u) = -2 7 + iu?J = D r {0) + d r (uj). (26) 

It would seem that since D r (0) oc the limit — > 00 
is problematic. However, as we will show, Z3 r (0) dis- 
appears from the equations for the spin Green func- 
tions. Since the bath is in equilibrium it must obey 
the fluctuation-dissipation theorem, so that D s (lu) = 
2icoth(cj/2T)ImL) r (w). Thus 

£ s H=-4 l7 ^coth(^). (27) 

This result can also be obtained from (|23(l with the 
Ohmic density of states, directly. 

IV. STATIONARY REGIME 

In this section we treat the retarded Green function 
G T (t,t') at large times t and t' , such that t, t' — ► 00, 
while t — t' is finite. We show that it is stationary, and 
that it is related to the stationary part of the statistical 
Green function via the fluctuation-dissipation theorem. 

The Lagrange multiplier z(t) is simply a function of 
one time variable, therefore for large t it should saturate 
to a constant value z^ ~ lim^oo z(t). We assume that 
in the coarsened phase z^ is independent of tempera- 
ture, which indeed turns out to be consistent with our 
findings. Since D T is independent of temperature, from 
Eq. ljTU|) it is evident that once z(t) saturates to z^, 
G T (t,t') no longer depends on temperature. This agrees 



with the classical result found in 0. Since G r {t,t') does 
not depend on G s , which is expected to age, one can 
find a time-translational solution for G r , once z x is ap- 
proached. This is a peculiarity of the p = 2 model, with 
the linear spin-bath coupling, since the retarded Green 
function is independent of the statistical Green function. 
We note that, consistent with our present discussion, no 
aging was found in the response of the infinite-range M- 
component rotors in the limit M — > 00 [5(- 

The time-translational invariance of G r for large times 
can be understood from a heuristic point of view. Be- 
cause the response function is decoupled from its statis- 
tical counterpart, it is independent of T, and therefore 
must behave in the same manner for all temperatures. 
At very large T, the system is paramagnetic and thus 
the response is time-translationally invariant. Thus it 
must exhibit this behavior at all temperatures. By con- 
trast, the statistical Green function, that is proportional 
to the spin correlation, is temperature-dependent so that 
the above argument does not apply to it. 

Assuming large enough times so that z(t) has satu- 
rated, and therefore G r (i, t') — G r (t — t'), one can Fourier 
transform Eq. (|10() over t — t' to obtain 



w 2 -z-d r ( w ) - J 2 G'-H' 

where z = z^ + Z3 r (0). This is a quadratic equation 
for G r (w) which also applies to the paramagnetic regime. 
However, z is still unknown, and will be obtained from 
the equation for the statistical Green function. It will be 
shown that z is finite, so that the divergence of D r (0) for 
51 — > 00 does not pose a problem. 

In order to obtain an equation for Gg T , we write the 
statistical self-energy in terms of the stationary and aging 
contributions, S s (i, t 1 ) = E| T (t - t 1 ) + T, s AG (t, t'), where 

£!t(*-*') = D s (t-t') + J 2 G s ST (t-t'), (29) 
Sag(M') = J 2 G s AG (t,t'). (30) 

Assuming t' is large, we substitute (|13|) into Eq. and 
obtain 



m-^ + Zoo ) [-2iq + G s ST (t - t')] 



At" [Eg T (t - t")G & {t" - t') + Y?{t - t")G s ST (t" - t')\ - 2iq 2J 2 G r (0) + £>■"(()) 

1 

I 



(31) 



In the derivation of Eq. H31[l. one encounters terms of the 
form J^df'G^it^'^Fit" - t'), where F can be any 

one of the stationary functions G r / a or lW a . We note 
that in these terms the aging function changes slowly 



compared to the stationary ones. We exploit this feature 
and assume that t" is close to t 1 , since the retarded and 
advanced functions decay to zero when their arguments 
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become large. We therefore approximate 

dt"G s AG (t, t")F{t" - t') ~ -2iqF(0), (32) 

where F(u>) is the Fourier transform of F(t). Similarly, 
dt"F(t - t")G s AG (t", t') ~ -2iqF(0). (33) 



In addition, in obtaining Eq. (|31(l . we used the fact that 
the imaginary part of a retarded Green function at zero 
frequency is zero. 

Although we assumed that t—t' is finite, we can Fourier 
transform Eq. I|31|l since it applies to any t — t' in the 
t 1 — > oo limit. The Fourier transform of Eq. I|31l) is 



jmw 2 — z — d r {uj) — J 2 G r (w) + G a (w) } G^ T (w) = 



D s (u)G a (uj) - 2iq5{uj) 2J 2 G r (0) + z 



(34) 



Since lim^oo G s ST (t) = the expression for G s st (lo) can- 
not include terms that are proportional to 8(u>). We ex- 
pect the other Green functions in Eq. (|34(l to be similarly 
well-behaved. Therefore in the coarsened phase, since q 
is finite, the term in the square brackets that multiplies 
the delta function must be zero, leading to the following 
equation for z and G r (0) 



-2J 2 G r (0). 



(35) 



In the paramagnetic phase q = 0, and therefore Eq. I|35|) 
is no longer valid. 



Another equation for z and G r (0) is found by substi- 
tuting uj = into Eq. fflfy. 



J 2 G r (0), 



G r (0) 



(36) 



which when compared with Eq. (|35|l renders [jG r (0)] 2 = 
1. Since G r must be analytic in the upper half plane, one 
obtains G r (0) = — 1/J and therefore 



z = 2 J, 



(37) 



in the coarsened phase. Indeed as we anticipated, z 
is a constant and does not depend on the tempera- 
ture T. Furthermore, z is finite, so the divergence of 
limn— »oo -D r (0) is not a problem; see the definition of z 
following Eq. and the discussion following The 
same z has been found for the classical p — 2 model ^} • 
In Ref. 2] an equation for z^ in terms of G r and £ r 
was found for the quantum disordered p-spin spherical 
models of p > 3; there E r depends on both G s and G r so 
that a simple solution for z^ does not emerge. However 
in the p = 2 case with linear spin-bath coupling E 1 is 
independent of G s . This decoupling leads to two distinct 
equations for Zoo and G r , and thus to solutions for both. 
More specifically (|35|l , the analogue of the equation that 
was found for in Ref. Q, contains only G r ; i(2"5|) pro- 
vides a second equation for z and G r , thereby permitting 
us to determine z. 

Now that z is known, we can solve Eq. (|28|) for G r (w) 
in the coarsened phase 



G r H 



2J 2 { 



mw 2 — 2 J + 2i"fu — [(mtij 2 + 2i^ui) (muj 2 

I 



4 J 



2H-)] 1/2 } 



(38) 



In solving the quadratic equation for G t (uj), we choose 
the negative sign in front of the square root to ensure 
that lim^^oo G r (oj) = 0. 

We now demonstrate the FDT in the coarsened phase 
for the stationary part of the Green functions. Using 
Eq. in Eq. , after the term with the delta func- 
tion has been eliminated, we obtain 



G s S tH - 



D s (lu) G t (lu) 



1 - J 2 G r (w) 



(39) 



We write D s (ui) in terms of InxD r (w) via the FDT, and 
taking the imaginary part of Eq. I|28|l we obtain 

G S ST H = 2 lC oth (^)lmG r H. (40) 



Thus, we have shown analytically that the stationary 
part of the statistical Green function is related to the 
imaginary part of the retarded Green function via the 
fluctuation-dissipation theorem, where the temperature 
is that of the bath. 

Now we have fully solved for the system's behavior in 
the stationary regime. Once the retarded Green function 
has been found, 138fl . the stationary part of the statistical 
Green function is known as well via the FDT, (I4UI) . 



V. THE EDWARDS-ANDERSON ORDER 
PARAMETER 

The Edwards-Anderson order parameter q can be 
found by substituting the expression that was obtained 
for the stationary part of the statistical Green function, 
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(|40|l (FDT), into the matching condition at zero time 
difference lfT7jl: 



^coth(^)lmG r H. (41) 



The Edwards- Anderson parameter is in the range < 
q < 1. We expect to obtain q = at equilibrium for 
the system in its paramagnetic phase. We note that 141(1 
admits negative values of q, which correspond to an in- 
stability of the coarsened phase; in this case one should 
look for a paramagnetic solution. 

As a check, expression (|41|) for q can be evaluated in 
two known limits. The first is at high temperatures, when 
T is much larger than all energy-scales apart from J. In 
this case coth(w/2T) ~ 2T/u and we obtain 

rji p DC J rp 

q = 1 + - / — ImG r H = 1 + TReG r (0) = 1 - -, 

7T J_ oc UJ J 

(42) 

where we have used the Kramers-Kro nig relation. This 
is the well-known classical result ^3 E3| > which should 
indeed be obtained when T is large. As an aside, we 
note that T should be smaller than J, since otherwise 
the system would undergo a phase transition into the 
paramagnetic phase. 

The other limit that we can check is T — > 0, since 
limT->o coth (lv/2T) = sign(w); in this case 



9 = 1 



^signHlmG». (43) 



In equilibrium, dw/(27r)sign(w)ImG r (w) = —1 |27j 
and we obtain q = 0, as expected. 

In our calculation of q we will use the following ap- 
proximation 



coth (—) 
\2TJ 



sign (a;) for 2T < lo, 
2T/lu for 2T > uj 7 



(44) 



breaks down, and the system must be in the paramag- 
netic phase. We have checked numerically in the 7^0 
and the T = planes, that there is no coexistence of the 
paramagnetic and the coarsened phases, except on the 
q = line. Furthermore, in section IVI Bl we show that 
the paramagnetic solution changes continuously into the 
coarsened one. Thus this is a second order phase tran- 
sition, and the q = line that is found in the coarsened 
phase is indeed the phase boundary. 

In the next section we describe the phase diagram. 
Here we concentrate on finding the quantum critical point 
which resides in the T = plane at the 7 — > limit. We 
expect the system to behave in a coarsened manner for 
larger values of to. Upon decreasing m we anticipate 
that the quantum fluctuations will grow until the mass 
reaches a critical value to*, where the system undergoes 
a phase transition into the paramagnetic state. 

At this point we comment on the physical meaning 
of taking 7^0. If one were to take 7 = from the 
start, one would initially decouple the spin system from 
the bath. In this case the FDT relationship that we had 
found, l|40|) . that includes the bath temperature, would 
not be valid — the spin system would not be "aware" of 
the bath's temperature unless it was coupled to it at some 
point in time. By taking 7 — > we consider the spin 
system to be coupled to the bath, with an infinitesimally 
small coupling. 

Substituting 7^0 into Eq. (J3BJ, we find 



ImG r (w) 



7-»0 



1 / TO \ 
J \j) 



1/2 



rriLU 
~4J 



1/2 



for 



M 2 ( — 

. TO 



1/2 



(46) 



and otherwise. Substituting |@SJ into Eq. for T 
0, one can calculate q exactly 



so that 



9 = 1- 



(47) 



2T 



2T 



q~l + — / — ImG r (w)+- / dwImG r (w). (45) 

7T J LU 7T J 2T 

We are aware that the model we are using for a bath can 
lead to peculiar results since it does not incorporate an 
infrared cutoff for the bath frequencies. In such a case 
the first integral in Eq. 14511 might be finite even when 
T — > 0. We have verified that this problem does not 
arise in our calculations, since for to — > ImG r (w) cx w Q , 
where a > 1. 

For a given average spin-spin interaction J we look 
at the phase space that is defined by the parameters 
to, 7, and T. We focus on the 7 — > 0, and the T = 
planes. The coarsened phase is defined by a finite q, 
while q = at any temperature corresponds to the para- 
magnetic phase. We shall look for the q = line, know- 
ing that beyond it (when q < 0) the coarsened phase 



Taking q = in l|47l) and solving for m we find the quan- 
tum critical point 



(48) 



We observe that this quantum critical point is identi- 
cal to that found for the model of M — > 00 infinite-range 
quantum rotors [l2T.ll3j]. This was originally identified us- 
ing replica methods which are tricky for quantum phase 
transitions since the ratio of the replica index to the tem- 
perature (M/T) must be kept fixed while M — ► 00 and 
T — > 0. It was also found using an asymptotic dynamical 
treatment, and it is reassuring that we recover the same 
result in the limit of zero dissipation for a treatment of 
a more general Hamiltonian. 
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IPM 



CPM 




CC 



QPM; 



T 

QC 



■-'(■'— 
J" 1 



subsections IVI Al and IVI 51 include our calculations, they 
are rather technical. For a summary of our results, we 
refer the reader to Tables [I] and [n] which are depicted 
graphically in Figs. [21 El an d E] The main qualitative 
features of the system's phase behavior are summarized 
in section IViTI 



A. The coarsened phase 

In the coarsened phase z — 2 J, and q varies; In Table [I] 
we summarize the expressions obtained for q in the 7 — > 
and T = planes. We identify the q = line as a second 
order phase boundary, where the paramagnetic solution 
changes continuously into the coarsened one, as it will be 
shown in section IYl Bl 



FIG. 2: The phase diagram in the 7 — + plane. The solid 
line depicts the second order phase transition line between 
the paramagnetic and the coarsened phases. One dashed line 
depicts the border between regions that are dominated by 
classical dynamics (CPM and CC), and the intermediate ones 
(IPM and IC) that are governed by both classical and quan- 
tum fluctuations. The other dashed line represents the bound- 
ary between the intermediate paramagnetic and the quantum 
paramagnetic (QPM) regimes, where the latter is dominated 
by quantum dynamics. The quantum coarsened (QC) regime 
exists only at T = 0, for m > m* , the dashed-dotted line. 



VI. THE PHASE DIAGRAM 

Using q and z as probes we now characterize the phase 
space that is spanned by the variables 7, m and T, by 
identifying the coarsened and the paramagnetic regions. 
When the system is in the coarsened phase, z — 23 [see 
(|o7fl ], and q ^ varies. The q = line in the 7, m, T 
phase space tells us when the coarsened phase becomes 
unstable. In the paramagnetic phase q — since the 
system is in equilibrium, and z varies. The spin-spin in- 
teraction is important in the coarsened phase, and indeed 
z depends on J. In the paramagnetic phase we expect z 
to depend on the two important energies 1/m and T. 

The system's phase behavior is presented in two sec- 
tions, [^^] and ^^Bj associated with the coarsened and 
the paramagnetic phases respectively. In the coarsened 
phase z = 2 J, and we characterize the behavior of q. By 
contrast, the system is in equilibrium in the paramag- 
netic phase, so that q = and we look for z. Ultimately 
we show that z changes continuously from the paramag- 
netic value to the coarsened one. The phase behavior in 
the 7-^0 plane is displayed in Fig.|3 and in Fig. |3we 
present our findings in the T — plane. In Fig. 01 we 
show the transition lines in the T — to plane for several 
values of 7, determined by a numerical calculation. Since 



1. The 7^0 plane 

Results (A) and (B) in Tableware obtained for 7 — > 0, 

by substituting the exact expression l|46[) for ImG r , 

7=0 

into Eq. 145|) . and performing the integration. This part 
of the phase diagram is depicted in Fig. |2 Expression 
(A) obtained for T > yj J/m is the well-known classical 
result, H42[) [T7L ITsj . which depends only on the tem- 
perature T. When T — and there are only quantum 
fluctuations, one obtains expression (C) in Table [I] for q 
[see also (03], which depends only on m. Expression (B) 
that was obtained for T < yj J/m has two competing 
terms: one that is governed by classical dynamics and 
resembles the classical expression (B), and another that 
is dominated by quantum fluctuations and is similar to 
(C). Therefore we call this regime an intermediate region, 
in which classical and quantum dynamics compete. 

We identify the line T c (m) that corresponds to q = 
in the m — T plane, as the phase boundary between the 
coarsened and the paramagnetic phases (c.f. solid line 
in Fig. [2J. For rn > 1/J the phase transition line is 
given by T c — J, independent of m. For m < 1/J it 
is found by setting (B) equal to 0. The line T c (m) is 
continuous approaching the point (to = 1/J, T = J) 
from both directions. Its derivative dT c /dm is contin- 
uous at this point as well, where it becomes zero. For 
to < 1/J the transition line terminates at the quantum 
critical point (m = m* , T = 0) with an infinite derivative 
dT c /dr7i — * 00. 



2. The T — plane 

The T = plane, depicted in Fig. |21 is by definition 
entirely quantum. In the weak coupling limit, 7 2 <C to J, 
we approximate 
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TABLE I: The Edwards- Anderson parameter q in the coarsened phase, for different regions of the phase space spanned by 
m, 7 and T; recall that here z = 2 J. The region labels CC, IC, and QC are according to Fig. |2]and Fig.|3] 



(A) 



(B) 



(C) 



(D) 



(E) 



m, 7 and T 



7^0 



7^0 

T < y/J/m 



7^0 
T = 



T = 

7 2 < mJ 



T = 

7 2 > mJ 



g(m,7,T) 



T 
7 



CC 



1 - 



2\l/2 



37TV mJ 



1-a; 



where x : 



„ a; \ 2 T . _i 
1 + i ^J Sm X 



j 



3n\/mJ 



21 7 



16 



ij 32 mJ mJ 



1-^-1^-0(1 
Ittj mJ V 7 



IC 



QC 



' mJ 



QC 



QC 



■'■T?('-f)K 



for < lu < 



2l 



1 / moj 2 I ^ iiij:- 



ImG r (w) ~ < 



J V J 
1 

J 2 

2~ 



4J 



2 \mujJ 



1^ , 2 7 /2J 
for — <u/< W— , 
J z m V rn 



(toJ + 7 2 )^ 



(tow 2 ) 2 70; 
" 4 + J 2 ' 



/2J AJ 4 7 2 

for W — < u)< \ h — ^ 

V to V to to 



/4J 4 7 2 
for a/ h — T < W. 

TO TO 



(49) 



Substituting these expressions into Eq. gSJ) we obtain (D) in Table □ 

The to = case is the extreme limit of strong coupling, 7 2 3> raj . In order to demonstrate that there are no 
coarsened solutions with m = 0, we substitute the imaginary part of the exact expression for G r (w), (|38[) . with m = 
into Eq. gSJ. We find that 



r 1 1 7^ 
q = — lim in ► —00, 

w-*oo 27T7 J 

verifying that the system must be in the paramagnetic phase when m = 0. 
In the calculation of q for 7 2 3> to J and to ^ we approximate 

~^"p Im { v/2»7w(2rya; - 4J)| for < w < J/7, 

for J/7 < cj < 2B/m, 
for 2"f/m < cj, 



(50) 



ImG'M ~ < 

27W 
2 7 



(51) 



r 



from which we find (E) in Table [I] The dominant loga- boundary for 7 2 3> to J (c.f. solid line in Fig. [3J we set 
rithmic part agrees with (|50)l . In order to find the phase 
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QC 



T 

1 - 



0.8 - 




QPM 



m 



FIG. 4: The q — lines in the T versus m plane for 7 = 



FIG. 3: The phase diagram in the T = plane. The solid 
line depicts the second order phase transition line between 
the quantum paramagnetic phase (QPM), and the quantum 
coarsened phase (QC). 



q = in (E) and obtain 

m = -j-e- 2 ^. (52) 

Thus the phase boundary for large 7 corresponds to an 
exponentially small m. 

From Fig.|3it is clear that as the coupling to the bath 
7 is increased, m must decrease in order to cross the tran- 
sition line. The spin-bath coupling causes a decay of the 
quantum fluctuations, which in the absence of tempera- 
ture, are necessary for the system to make the transition 
from the coarsened to the paramagnetic phases. 



3. Finite T and 7 

In Fig. 0] we display the phase boundaries T c (m) for a 
set of spin-bath couplings 7. We determined these lines 
numerically, by solving q\ (m) — for T, at a set of m 

values, where q and G r are given by l(4"T|) and re- 
spectively. For constant m at low temperatures, as 7 is 
increased, T must be raised in order to destabilize the 
coarsened phase. As previously discussed, the coupling 
to the bath causes a decay in the quantum fluctuations. 
Therefore its effect can be overcome either by decreas- 
ing m, thereby increasing the quantum fluctuations, or 
by increasing T. As an aside we note that at high tem- 
peratures, increasing 7 has a much smaller effect. Our 
finding that the position of the transition line depends 
on the strength of the spin-bath coupling is qualitatively 
similar to that obtained in recent studies of dissipative 
effects on more complex quantum systems psf . 



TABLE II: The Lagrange multiplier z for different regions of 
the phase space spanned by m, 7 and T, in the paramagnetic 
phase. The region labels CPM and QPM are according to 
Figs. Hand El 





m, 7 and T 


2(771, 7, T) 


(F) 


7^0 


T 1 + (jpj CPM 


(G) 


7^0 


— [l + 12(mJ) 2 + O(mJ) 4 ! QPM 
Am 1 J 


(H) 


T = 

7 2 < mJ 


[l - — + 12(mJ) 2 + • • -1 QPM 
4m L 7r J 


(I) 


T = 

7 2 > mJ 


-^-exp(-27T7) QPM 
m 



B. The paramagnetic phase 

The expression for z that we found earlier, 1)37(1 . applies 
only to the coarsened phase. In order to characterize the 
paramagnet, we must find an expression for z in this 
regime. Therefore we consider Eq. I|45ll with G r given by 
the result of Eq. 1(28). substitute q = 0, and solve for z. 
In Table |n] we summarize the results of this calculation. 

We remind the reader that Eq. (|28|l for G r applies to 
both the paramagnetic and the coarsened phases, where 
their difference lies in their respective values of z. At the 
transition line, z = 2 J for both phases. We find that in 
the paramagnet z is larger than its coarsened value of 
2 J, and that it decreases continuously to its coarsened 
value at the transition line. For the cases of 7 — > 0, and 
of T = 0, we have checked numerically that an attempt 
to solve q = using ((2*S|l for G 1 ', could be obtained be- 
yond the transition line (in the coarsened region) only 
if z continued to decrease below the value of 2J . Such 
a solution would incorporate a finite ImG r (w) at uj = 
which is unphysical. Thus we deduce that the coarsened 
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and the paramagnetic phases coexist only at a single line, 
corresponding to a second order phase transition. 

Let us begin by finding z for large T. From Eq. Ij28(l 

"4/2)1 /2J 2 . For T 



we know that G r (0) 



z + VP 



much larger than all other energy-scales, including J, and 
for q = 0, we find that g = = 1 + TReG r (0) [following 
Eq. d2M ) and therefore 



T 1 



y2 



(53) 



in agreement with the classical result ^jj- Note that 
when T = J and the g = line is approached, we obtain 
z = 2J from J53J); e.g., for the specific case of 7 — ► 0, see 
the solid line for m > 1/ J in Fig. [21 Thus z, and therefore 
also G r (w) and G s (w) change continuously when q = 0, 
and indeed T = J is a second order phase transition line. 



1. The 7^0 plane 

Substituting 7 — > in Eq. (|28|) and solving for G r , we 
find that 



ImG r H = -jJ 1 



mw — z 



2 J 




for uj- < uj < cj+, 

. (54) 

where cj± = y/(z± 2J)/m, and ImG r = otherwise. 

If 2T > u+, we use (I45|l to calculate q, and then, sub- 
stituting q = we obtain the following equation for z, 



(55) 



The exact result of the integration on the left hand side 
is a combination of elliptic functions, from which it is 
hard to isolate z. We will therefore expand in 2 J/z. It is 
convenient to make a change of variables sin 6 — (mui 2 — 
z)/2J . This procedure leads to expression (F) in Table ITT1 
for z [see also l|53|)] . Substituting this expression for z into 
uj + we find the line T = lo + /2, the boundary between the 
IPM and the CPM regions above which the system seems 
to be predominantly classical (c.f. dashed line in Fig.[3J). 
For 2T < uj- , we obtain in the same manner as above 




If one neglects 2 J/z completely, the result is 

1 

4m' 



(56) 



(57) 



in agreement with This result is valid for m <C 1/8 J. 
In this case 2J/z <C 1 and our approximation is consis- 
tent. 

If we continue the expansion up to linear order in 2 J/z 
we obtain expression (G) in Table [HI Substituting m* 



into (G) we find z = 1.93 J which is very close to the 
expected value of z = 2J in the coarsened phase. Thus, 
it would seem that even close to the transition between 
the coarsened and the paramagnetic states this is a good 
approximation. Substituting (G) into the expression for 
uj- we were able to draw the line T — cj_/2, the dashed 
line between the QPM and the IPM regions in Fig. |3 
below which the Lagrange multiplier z is governed by 
quantum behavior. 

When cj_ < 2T < ui + the system is in an intermediate 
regime that is governed by both classical and quantum 
dynamics (IPM regime in Fig. We define 6t such 
that sin Q T = (AmT 2 - z)/2J. Using Eq. in order to 
calculate q and then setting q = we obtain 



2T 



dO- 



cos 2 9 



-n/2 



z + 2 J sin 
1 



(58) 



ny/m J 0T ■sfz + 2 J sin 9 
Upon expansion in 2 J/z, the equation for z becomes 



|4cos 3 T + (59) 

37TZ 




J 

37TZ 



cos 3 e T } = 1. 



Clearly the first term in Eq. I|59|l describes the classical 
fluctuations due to temperature, while the second de- 
scribes the quantum fluctuations due to the mass. 

This point can be further clarified by rewriting Eq. dS 
in the following form 



2T 



h 



mT 2 



1 (mT 2 z 
Imz \ J J 



1, (60) 



where /1 , ji are slowly changing functions of their vari- 
ables. Neglecting terms of the order of (J/z) 2 , for 
T = UJ+/2, we obtain f\ ~ 1/2 and fi — 0, while for 
T = u)-/2 we get /1 = and f 2 ~ 1/2. Solving Eq. ^ 
for z as if f± and f% are known we find that 



2m 



2/iT 



h 



fiT. 



(61) 



When T = uj + /2, Eq. JSTJl leads to z = T, while for T = 
uj-/2 one obtains z = 1 /Am; the results for T > lo+/2 
and T < ui-/2 up to the leading order, respectively [c.f. 
(F) and (G)]. 



2. The T = plane 

For T = it is convenient to begin our study in the 
paramagnetic region, far away from the phase boundary. 
Since the spin-spin interaction is not necessary for the 
system to be in its paramagnetic state, we can assume 
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at first that J 2 <C min {47 2 z/m, z 2 }, and neglect J in 
Eq. i|28|) . In this case the retarded Green function is given 
by G T (°>(u)) = (mui 2 + 2ijuj — z)~ l , where the superscript 
stands for zeroth order in J. Thus we obtain from 
Eq. (|43[) at q = 0, T — and J = the following equation 
for z 



{■y — \l mz — 7 2 
in — = 1 . 



I-Kisjmz — 7 2 «7 + y/mz — j 2 



(62) 



In the weak coupling limit, expanding Eq. I|62|) for 7 2 <C 
mz, we find 



1 

4m 



1 - 



8- 



(63) 



which for 7 — > is consistent with Eq. 

Now we reintroduce J iteratively, by considering 
G r W = (muj 2 + 2ijtv — z — J 2 ^ )) -1 . We obtain from 
Eq. (|43|) at g = and T = the following equation for z 



1 — - arctan 



y/m(z+J) — i- 2 



\Jm{z + J) — 7 2 

1 — — arctan 



(64) 



V m(z-J)-7 2 



\Jra{z — J) — 7 2 



1. 



Expanding Eq. I|64|) in 7 2 /mz and in J/z, since we expect 
J to introduce a small correction to z in l|63|l , we find (H) 
in Table ITU in agreement with (G) for 7 — > 0. 

Returning to Eq. I|62|l and assuming strong coupling 
so that 7 2 ^> mz, we obtain (I), which will equal 2 J 
when m — 2^ 2 / Jexp (— 2777). This is exactly the same 
expression for the transition line that we obtained for the 
coarsened phase [see (|52*)l ]. Thus for large 7 the transition 
is also continuous. 



VII. DISCUSSION 

In this paper we present a study of the phase behav- 
ior of the disordered quantum p = 2 spherical model 
as a function of temperature, dissipation (spin-bath cou- 
pling) and strength of quantum fluctuations (mass) at 
long times, identifying different paramagnetic and coars- 
ened phase regions. Using the weak-ergodicity breaking 
ansatz and a linear spin-bath coupling, we characterize 
the system's behavior in the stationary regime. We do 
so by exploiting the fact that its retarded Green func- 
tion is independent of its statistical counterpart, and is 
time-translationally invariant at long-time scales. In do- 
ing so, we demonstrate analytically that in the stationary 
regime the statistical and the retarded Green functions 
are linked by the fluctuation-dissipation theorem (FDT). 

We use the spherical constraint and the fluctuation- 
dissipation theorem in the stationary regime to charac- 
terize the phase behavior of the system. More specifi- 
cally, the long-time behavior of z, the Lagrange multi- 
plier that enforces the spherical constraint, is crucial for 



distinguishing different paramagnetic regimes. Similarly 
the Edwards- Anderson order parameter, q, determined 
from the stationary statistical Green function (FDT) 
and matching conditions at equal-times (spherical con- 
straint), characterizes different coarsened phases. We 
find that all the boundaries between the paramagnetic 
and the coarsened regions, determined by the q — lines, 
in the 7 — > and the T — planes are continuous; this 
is in marked contrast to the situation for quantum p > 3 
spherical models where quantum fluctuations at T = 
drive the transition first-order 0,H)@- Furthermore, we 
identify a quantum critical point for T = and 7 — > 
that coincides with the one found for a model of infinite- 
range, infinite-component quantum rotors in the absence 
of dissipation |3 . 

The Edwards-Anderson order parameter is used to 
probe the different types of coarsened phases; here z is 
a constant (z = 2J). More specifically, the behavior of 
q determines whether the system's dynamics are domi- 
nated by classical or by quantum energy-scales. From our 
results, it is clear that the increase of either quantum 
or thermal fluctuations (c.f. Figs. |2] and 01, by raising 
the inverse-mass or the temperature respectively, leads 
to the destabilization of the coarsened phase. Further- 
more an increase of the spin-bath coupling, 7, causes a 
decay of the quantum fluctuations; the inverse-mass must 
be enhanced still further to approach the phase bound- 
ary. Finally for constant mass and increasing 7 at low 
temperatures, an increase in T is necessary to cross the 
q = phase boundary (c.f. Fig. [3J|; the effect of the 
spin-bath coupling is significantly reduced at higher T. 

In the paramagnetic phase, q = 0, and the system 
is characterized by its long-time Lagrange multiplier, z, 
which enforces the spherical constraint. In both the 
7^0 and the T = planes (c.f. Figs. H and re- 
spectively) a quantum paramagnetic phase exists in the 
limit of vanishing mass. However in the 7^0 case there 
are crossovers between different paramagnetic regions as 
both temperature and mass are increased, corresponding 
to states where classical energy-scales assume increasing 
importance compared to their quantum counterparts. 

The simplicity of this quantum p = 2 spherical model 
is both its blessing and its curse. It is exactly this fea- 
ture that permits us to show explicitly the validity of the 
fluctuation-dissipation theorem in the stationary regime; 
furthermore this system is accessible to analytic study in 
the limit of vanishing dissipation. However, we should 
also point out that here the response is always indepen- 
dent of temperature, a feature that makes it particularly 
awkward for comparison with experiment in the classical 
and quantum regimes. Furthermore we do not have any 
spatial information, so possible overlap with real coars- 
ening measurements is rather difficult. However some di- 
alogue between theory and experiment may be possible 
near the quantum critical point, where issues of length- 
scales assume reduced importance. As a first step to- 
wards making such contact with experiment, we need to 
establish whether the quantum critical point studied here 
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is in the same universality class as that of the Ising model 
in a transverse magnetic field. It is intuitive to relate the 
latter to the M = 1 quantum rotor case 0, HD, Ifllj . 
which would suggest that it is in a different universal- 
ity class from the p — 2 quantum spherical model (that 
strongly resembles the M — > oo case) studied here. How- 
ever the distinction between continuous and discrete spin 
symmetries may be important |8(1 . so this argument is 
inconclusive. On another note, it is known that one can 
glean spatial information from the study of the classical 
p = 2 spherical model by exploiting a mapping on long 
time-scales to the three dimensional 0(n) model [l8Ll29| . 
and an analogous mapping should be explored for the 
quantum system. Finally we plan to characterize the dy- 
namics of this simple system near its quantum critical 
point in order to study whether it has any slow relax- 
ational modes that are not shared by its classical coun- 



terpart. 

Note added in proof: Recently, we learned of earlier 
studies [32l ] of the critical behavior of quantum rotors at 
finite dissipation. In the relevant parameter regime our 
results are qualitatively similar. We thank S. Sachdcv 
and A. Sengupta for bringing these papers to our atten- 
tion. 
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